Self-assembled aggregates in the gravitational field: 
growth and nematic order 



m 
O 
O 

(N 
Oh' 

<3\ 



Vladimir A. Baulin* 

Service des Interfaces et des Materiaux Moleculaires et Macromoleculaires, 
DRFMC, CEA-Grenoble, 17 rue des Martyrs, 38054 Grenoble, Cedex 9, France 

The influence of the gravitational field on the reversible process of assembly and disassembly of 
linear aggregates is focus of this paper. Even the earth gravitational field can affect the equilib- 
rium properties of heavy biological aggregates such as microtubules or actin filaments. The gravity 
gives rise to the concentration gradient which results in the distribution of aggregates of different 
lengths with height. Strong enough gravitational field induces the overall growth of the aggregates. 
The gravitational field facilitates the isotropic to nematic phase transition reflecting in a broader 
transition region. Coexisting phases have notedly different length distributions and the phase tran- 
sition represent the interplay between the growth in the isotropic phase and the precipitation into 
nematic phase. The fields in an ultracentrifuge can only reinforce the effect of gravity, so the present 
description can be applied to a wider range of systems. 
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I. INTRODUCTION 



Living polymerization has a substantial number of in- 
dustrial applications in the production of soaps, dyes, 
organometallic complexes, etci£ Moreover, some bio- 
logical objects owe their functionality to reversible pro- 
cesses of assembly and disassembly. Particular examples 
are microtubules, actin and intermediate filaments com- 
posing the cytoskeleton of living cellsi 3 ^ 4 ^ They are in- 
volved in several vital processes such as transport, chro- 
mosomal segregation,- 6 growth and division of a cell. In 
contrast to common synthetic living polymers, biologi- 
cal self-assembled objects have a considerable molecular 
masses, sufficient to induce a significant concentration 
gradient in an aqueous suspension of these objects. In 
turn, the concentration gradient can influence the equi- 
librium size of the objects and for long rigid aggregates 
it can even provoke the isotropic - nematic phase transi- 
tion. 

Before we address such phenomena it is useful to get 
an estimate of the relevance of the effect of gravity for 
a real system. Let us focus on a particular example of 
microtubules. The typical length of a microtubule is 10 
jj,mfL^ the linear density is 1.6 x 10 5 (g/mol)/nm, so the 
total mass of a microtubule is ~ 10 9 g/mol. The molec- 
ular mass of a tubulin dimer, subunit of a microtubule, 
is mo ~ 10 5 g/mol^ thus, one microtubule comprises of 
N ~ 10 4 tubulins. The influence of the gravity is signifi- 
cant if the potential energy of the object in the gravita- 
tional field is of order or higher the thermal energy, kT, 
i.e. if (mgh/kT)N = /3N > 1, where m — mo — pv is the 
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mass of a subunit corrected for buoyancy (jo ~ 1 g/cm 3 is 
the density of water, v ~ 100 nm 3 is the subunit volume, 
so the effective mass in a solution is m ~ 10 4 g/mol). 
Thus, the gravity is important when [3 > /3* = 1/N. 
Assuming the total height of the vessel, h ~ 1 cm, the 
gravitational acceleration, g ~ 10 m/s 2 , the temperature, 
T ~ 300 K, we get /3 ~ 10/3*. So, even the earth gravi- 
tational field can influence the equilibrium properties of 
such heavy objects as microtubules. Note, that one can 
reach 10 4 — 10 5 times larger gravitational acceleration in 
an ultracentrifuge, 9 which sometimes is used in the pro- 
cess of preparation of microtubules, 8 or in the pelleting 
experimentsii2*ii To this end, the gravitational field can 
serve as a driving force for the assembly and, eventually, 
to the nematic ordering of growing aggregates. Another 
way to study the equilibrium distribution of linear aggre- 
gates is the rotating clinostat experiments 12 which con- 
stitute the rotation of a tube about the horizontal axis. 
Changing the direction of the "fall" all the time can ef- 
fectively reduce the gravity. Such apparatus gives the 
gravitational fields in the range 10 -2 -j- 10 2 g. 

The rotating clinostat experiments 1 ^ were used to ex- 
amine the microtubule self-organization and structura- 
tion by reaction-diffusion mechanisms leading to different 
pattern formation. It was shown, that microtubules self- 
organize in different patterns in the normal gravitational 
field, while no self-organization occurs in the absence of 
the gravityi 13 i 14 i 15 

It is noteworthy, that the gravity can also bring about 
the opposite process, namely the floating to the surface 
when m < pv. One of the most known processes of that 
type is the creaming of emulsions. In our consideration 
this only changes the sign of the effective mass, m, while 
the top and the bottom switch their places. 

The problem of the influence of weak gravitational 
fields on different vital processes is largely investigated 
in the framework of so-called microgravity research per- 
formed in sounding rockets, parabolic flights or a space 
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station. In particular, it concerns a number of experi- 
ments conducted on board of space-shuttle flights. For 
example, the formation of liposomes is altered by the 
gravitational fields the size of the liposomes formed in 
the absence of the gravity is considerably larger than 
that formed in the normal earth conditions. The gravity 
also influence the long-range ordering of some systems: 17 
space-grown crystals have better optical quality, smooth- 
ness of crystal faces, etc. The problem of crystal growth 
under gravity is also discussed in ref. llSl 

However, the number of theoretical treatment of these 
observations is not so impressive. Here are some exam- 
ples: the size distribution with height of noninteracting 
linear micelles in the gravitational field has been analyzed 
in ref. Il9i It was shown, that the concentration gradient 
due to gravitational or centrifugal field leads to the hy- 
perbolic dependence of the average aggregation number 
with height in contrast to simple barometric dependence 
for nonaggregating objects. The dynamics of linear mi- 
celles in the ultracentrifuge is studied in ref. 20. The 
sedimentation of colloidal rods in an ultracentrifuge, has 
been studied theoretically in ref. l2ll 

In our previous papen 2 ^ we studied the isotropic to ne- 
matic phase transition of rigid monodisperse rods with 
fixed size in the gravitational field. The major result is 
that the gravity can induce the nematic ordering in the 
system of rigid rods, while the main effect of the grav- 
ity on the phase diagram is the broadening of the phase 
coexistence region. It was also shown that interactions 
between rigid rods can lead to deviations from the baro- 
metric distribution of rods with height. 

In the present paper we analyze a wider range of sys- 
tems. Instead of fixed length of rods^ we assume that 
the size of rods determines by the equilibrium conditions. 
In particular, this is the case of (i) living polymers, which 
can attach monomers to the end of the chain, (ii) linear 
surfactant micelles, (iii) biological self-assembled objects 
such as microtubules and actin filaments. Although all 
these systems have different geometries and different size 
scales, they have additional degree of freedom, namely 
the annealed length of the aggregates which can be tuned 
by external conditions. Thus, we expect, that the gravi- 
tational field can have a severe influence on the equilib- 
rium size and the concentration profile of the aggregates. 
The new result of this paper is that the gravity induces 
the overall growth of the aggregates. 

Concerning the isotropic to nematic phase transition of 
stiff aggregates with annealed size, we expect the same 
effect of the gravity as in the case of rods with fixed 
length. Namely, the gravity can induce the nematic or- 
dering and broaden the coexistence region. However, an- 
nealed length implies a large polydispersity of the ag- 
gregates, thus, we investigate the influence of the grav- 
ity on the distribution of lengths in a solution. Since 
the phase separation is coupled with the redistribution 
of aggregates of different lengths with different heights, 
we discuss the structure of the two coexisting phases and 
the distribution of the average aggregation number with 



height. 

The paper is organized as follows. In the next section 
the general definitions and the expressions for the free 
energy and the distribution function are given. In the 
section III we consider the equilibrium properties of the 
homogeneous isotropic solution of living polymers under 
gravity. The concentration profiles and the average ag- 
gregation numbers as a functions of height are analyzed. 
In the section IV, we examine the isotropic to nematic 
phase transition in the gravitational field. We discuss 
the effect of the broadness of the phase separation region 
in the gravitational field for living polymers, the length 
distribution functions of the aggregates in the coexisting 
region and polydispersity of the phases. In the Conclu- 
sion we summarize the obtained results. 



II. FREE ENERGY AND DISTRIBUTION 
FUNCTION 

Living polymers are linear aggregates which can bind 
reversibly one subunit to the end of the chain. They 
show a broad distribution of lengths, since their energy 
depends only on end effects and there are no packing 
constrains preventing the infinite growth as in the case 
of spherical micelles. The scission end energy (in units of 
kT), 5 > ensures energetically favorable formation of 
aggregates. The aggregation energy per subunit, ejv has 
a simple formi^ 3 . en = £oo + S/N, where is the bulk 
energy, which is independent on TV. Thus, Nsn is the 
energy of an aggregate with the aggregation number TV. 
The subunit concentration is assumed to be much larger 
than CMC, so the aggregates are rather long. 

Our description is based on the theories of living poly- 
mers and the liquid crystalline transition in the absence 
of the fieldi 24 . 25 i 26 i 27 i 28 . 29 i 30 i 31 i 32 i 33 i 34 For our descrip- 
tion we take as a starting point ref. l30l It is known, 
that the assumption of perfectly rigid aggregates leads 
to unphysically strong coupling of growth and nematic 
ordering. 30 . 32 . 35 

That is why, we have to allow for the partial flexibility 
of the aggregates. We do it in accordance to refs. 32,33 ' 35 , 
where the flexibility of the chains is coming through a 
particular form of the trial function in the nematic phase. 
However, in this approach the flexibility of the isotropic 
phase is not taken into account. 

We consider the solution of living polymers in the ves- 
sel subject to the gravitational field. Suppose, that the 
gravitational field acts in the z direction, such that the 
subunit concentration varies from its maximal value at 
z = (bottom) to its minimal value at z = 1 (top). Our 
discussion is focused on the case of long aggregates. We 
ignore the end effects, though they can be important for 
short chains^ , that dominate at the top of the vessel, but 
for the subunit concentration well above the CMC, the 
fraction of short aggregates in a solution is low. 

In the following we will see, that the gravity can lead 
to the nematic ordering in the described system, thus, 
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we write the free energy in a general form which assumes 
the orientational order. The free energy can be written 
in terms of the dimensionless concentration of rods each 
comprised of N subunits oriented in the spatial angle f2 
on the height z as C(N, fi, z) (this is the number concen- 
tration multiplied by the volume of a single subunit, v). 
C(N, fi, z) obeys the following constraint 



NC(N, n, z)dNdVL = <j>{z) 



(1) 



where 4>{z) is the total volume fraction of subunits on 
the height z. The total subunit concentration is $ = 
J (j>(z)dz. 

The free energy is composed of an ideal and an inter- 
action terms, F = Fid + Fi nt . The ideal term per unit 
volume is written in the form 



F ld {z) 
kT 



C{N, n, z) [In (4irC(N, ft, z)) - 1+ 



Ne N ]dNdn + J C(N,n,z)(3NzdNdn(2) 

The first term is the well known expression for the free 
energy of self-assembled aggregates^ The factor Air is 
introduced for convenience of description of the isotropic 
phase, though it can be absorbed in the energy of an ag- 
gregate, jVejv- The last term is the potential energy in 
units kT of aggregates in the gravitational field. The only 
parameter associated with the field is (3 = (mgh)/kT, 
which depends on the total height of the vessel. Note, 
that although the last term in (J2J coincide with that in 
the case of monodisperse rods of fixed length^ the defi- 
nitions of this parameter does not coincide. Here m is the 
mass of a single subunit, while in ref. |22| m is the total 
mass of a rod. Again, we assume that the characteristic 
length of the variation of the field is much bigger than the 
length of a rod, although the fields in the ultracentrifuge 
can be very high and this assumption can be violated. 

In order to apply the Onsager approximation^ 3 ^^ we 
assume excluded volume interactions only. Moreover, 
only pair interactions are taken into account, higher or- 
der interaction terms are neglected. This is a strong ap- 
proximation which is strictly speaking not valid for high 
concentrations in the nematic phase. As we will see, the 
strong gravitational field induces the precipitation of ex- 
tremely long, and thus, heavy, aggregates and the second 
virial approximation should be corrected. 

Nevertheless, this approximation works well in the 
isotropic phase and it also gives reasonable qualitative 
estimates in the nematic phase. 

The interaction term per unit volume for athermal ag- 
gregates, written in the second virial approximation has 
the following form^S^ 



Fir ^ = - [ C(N,n,z)C(N',fl',z)B 2 h)dNdN'dndn' 
kT v J 

(3) 



where -82(7) is the second virial coefficient which de- 
pends on the angle 7 between two rods. Neglecting the 
end effects, it can be written in the fbrm25*2& -62(7) = 
LL'D |sin 7I , where L and L' are the lengths of two cylin- 
ders, D is their diameter. We can rewrite this expression 
in terms of aggregation numbers N, N 1 and the volume 
of a subunit, v, as -62(7) = xvNN' |sin'y|, where x is 
the parameter associated with the form of a subunit. In 
the case of cubic subunits, h = 1, while for cylindrical 
subunits, x = 4:1 /D, where I is the height of a cylinder. 
However, to apply this notation to the case of micro- 
tubules, we should note, that a microtubule has complex 
geometry. Microtubules are hollow cylinders consisting 
of / = 13 columns of tubulin on the surfaced So, they 
grow up at / places at the same time. Thus, the length of 
a microtubule L ~ fN and x — Alf 2 /D. In the following 
we assume k= \ for simplicity. 

Variation of the free energy J2J), © with respect to 
the concentration C(iV,f2, z) along with the constraint 
condition leads to the following expression 



C(N,n,z) 



exp(JV(ju — en — flz)- 



1 

47T 

~ f C(N,n,z)B 2 (j)dNdn\ (4) 



where [i is the Lagrange multiplier, which can be elimi- 
nated with the help of the constraint JIJ. 

Introducing the average aggregation number of rods in 
the direction 17 at height z as 



N(fl,z) 



J NC(N, Q, z)dN 
J C(N, O, z)dN 



(5) 



we can write the expression for the distribution function 
C(N, fi, z) in a simpler form 



jN 2 (si,z)dn 



■ exp 



N 



N(n,z) 



(6) 



To write the last equality we used the constraint condi- 
tion (HJ). 

The final expressions could be considerably simpli- 
fied if we note, that (i) there is no angular depen- 
dence in the isotropic phase, (ii) in the nematic phase 
we can split the distribution function C(N, f2, z) in 
two parts^ one dependent on the aggregation num- 
ber, N, and another dependent on spatial angles, tt. 
Namely, we denote the length distribution function, 
P(N, z) = j NC(N, f2, z)d£l with the normalization con- 
dition J P(N, z)dN — cj)(z) and the orientation distribu- 
tion function R(N : Q, z) = j Nc(N'n 'z)dn w ^fi the nor- 
malization / R(N,n,z)dtt = 1. With this, C(N,n,z) = 



4 



^^-R{N, n, z) and the free energy, F = F id 
now yields in the form 



Pii 



F(z) 
kT 



P(N, z) 
N 



dN 



a(N,z) 



1 P (^ z ) i 



Ne N + l3Nz + -N \ P(N', z)p(N, N', z)dN' 



(7) 



where a(N,z) = f R(N, f2, z) In AttR{N, f2, z)d£l is the 
orientational entropy, p(N, N',z) = — f R(N 7 0, z) 
R(N', ft', z) |sin7| dCldil' is the excluded volume term. 
We keep this nomenclature in order to be consistent with 
refs i 22 ^? 2 ^?!?^ In the isotropic phase ai = 0, pi = 1 
and Ri(N,Cl, z) = 1/4tt, so the expression for the free 
energy simplifies significantly. 

In the next section we describe the concentration pro- 
file and the average length of the aggregates assuming 
that the solution of rods is homogeneous and isotropic 
for the whole range of z. 



III. HOMOGENEOUS SOLUTION: 
OF AGGREGATES 



GROWTH 



<N> 



1200 n 




1100 



1000 



0.000 



0.001 



0.002 



0.003 



0.0O4 



FIG. 1: Average aggregation number on the whole range of 
heights, (Ni), of ideal linear aggregates 11411 as a function of 
a dimensionless parameter associated with the gravitational 
field, P = mgh/kT (solid line) in comparison with the ex- 
act solution for nonideal aggregates (squares). The parame- 
ters used: 8 — 25, $ = 10~ 5 (low concentration limit). The 
crossover between the week field and the strong field, which 
induce considerable growth of aggregates is /3* ~ l/iVi.o. 



Since the directions of rods in the isotropic phase are 
completely disordered, the distribution function does not 
depend on angles, Cj(iV, Q) = Ci{N) and eqs. Q), © 
give^i the well known result for the average aggrega tion 
number in the isotropic phase Ni(z) — \J 4>i(z)e s . In 
turn, the free energy expressed in terms of the length 
distribution function, P(N, z) J7J) leads to the free energy 
of the isotropic phase in the form 



kT 



Pz- 



Ni{z) 



(8) 



Here we used again the phcnomcnological expression for 
the aggregation energy (per kT) of linear aggregates^ 
£jv = Soo + S/N. The chemical potential (ii(z) = 
dFi(z)/dMz) is 



all aggregates have almost the same potential energy. In 
our approximation the gravitational field influence the 
system only through the concentration gradient. The ex- 
pressions © , II 01) will be of use to the phase coexistence. 

In the equilibrium the chemical potential does not de- 
pend on height. Thus, p,i(z) = const, eq. determines 
the concentration profile induced by the field. In turn, 
the variation of the concentration with height leads to 
a redistribution of aggregates with different lengths at 
different heights. To address this issue, we focus on two 
limits: ideal aggregates (low concentrations) and consid- 
erably long aggregates (high concentrations). 

In the limit of noninteracting ideal aggregates in the 
gravitational field^ the last term in the expression for 
chemical potential © arising from the excluded volume 
interactions is eliminated. This leads to the following 
concentration profile 



kT 



■j3z- 



Ni{z) 



(9) 



and the osmotic pressure Pi(z) = 1/vc/r-d [i^/^i] /d<j>i{z) 
is 



= V -d>i(z) 



(10) 



Note, that these expressions formally have the same 
form as in the absence of gravity^ ^ (the only differ- 
ence is that the concentration 0j depends on z). This is 
because on the lengths smaller than the variation of grav- 
itational field, the system does not "feel" the field, since 



<f>i( z ) = 



I + y/<l>i iO e s 0z 



(11) 



where — 4>A Z — 0) is the concentration at the bot- 
tom. This expression is interesting because it shows non- 
barometric dependence as expected for the case of nonag- 
gregating objects ml 

The average aggregation number, Ni(z) = \J 4>i(z)e 5 , 
has the hyperbolic dependence on the height z*& De- 
noting the average aggregation number at the bottom, 
Ni, = Ni(z = 0),it yields 



Ni{z) 



N 



t.O 



1 + N ifi Qz 



(12) 
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FIG. 2: Typical concentration profile of living polymers sub- 
ject to the gravitational field in a dilute solution (thick line) 
for f3 — 0.7, S = 11 and the total concentration "!>i = 0.1. 
The dependence 4>i(z) is linear except the region on the top 
where the concentration is small. The crossover between the 
regimes is z cr — \J Qiiv/fl. In contrast, the concentration pro- 
file of nonaggregating subunits (S — 0) for the same total 
concentration and the field (thin line) obeys the barometric 
distribution. 



This expression shows, that the aggregates at the bot- 
tom are longer than that on the top. This is simply a 
consequence of the redistribution of subunits with height. 
It does not tell about the change of the equilibrium size 
of the aggregates. In order to examine the effect of grav- 
ity on the system as a whole, one can consider the global 
average on the whole range of heights, (Ni) , as a function 
of j3 = (mgh)/kT. This average is related to the total 
subunit concentration, $^ = J (f>i(z)dz, which assumes to 
be a constant 



(Ni) = 



J NC t (N,z)dNdz 
J d{N, z)dNdz ~~ jN z (z)dz 



(13) 



If the influence of the gravitational field is rather weak, 
i.e. when (3 <C (3* ~ it coincides with the expres- 

sion in the absence of the field, (Ni) ~ \J <&ie s . Thus, low 
gravitational field induces only the redistribution of ag- 
gregates with height and it does not influence the micel- 
lar equilibrium properties. On the other hand, if (3 > /?*, 
(Ni) yields^ 



(N) = 



hi 



(14) 



This is an increasing function of (3, i.e. strong enough 
gravitational field is able to induce considerable aggrega- 
tion even of ideal micelles. The typical plot of this func- 
tion is presented in Fig. 1. Note, that this expression 
is valid only for extremely low concentrations, <C 1. 



FIG. 3: The average aggregation number of living polymers in 
a dilute solution, Ni(z), as a function of the height z for the 
same parameters as in Fig. 2. The plot can be approximated 
by two asymptotes: long aggregates, interacting via the sec- 
ond virial approximation for z < z cr , and ideal aggregates for 

Z > Z cr with Zcr — \J QiTt/fl. 



However, the deviation of this expression from the exact 
solution (with the interaction term) is rather small. We 
can also calculate the asymptote for (3 <§; 1, 



{Ni) « + - (^e s f 2 /3 2 + 0(f3 4 ) (15) 

If the concentration is not low, we can still get the 
analytical result in the limit of long aggregates, i. e. when 
we can neglect the third term in eq. This can be 

done for high enough end energies 5, i.e. in the formal 
limit, i > 1. In this case, the concentration depends 
linearly on the height, 



(f>i{z) = <j> ii0 f3z 



(16) 



while the average aggregation number yields in the form 



Ni{z) 



N 2 _ Z e sp z 



2 

— < 

7T 



(17) 



where N^q = \J e (5 ($i + /3/tt) is its value at the bottom, 
z = 0. This expression brings us immediately to the con- 
clusion, that the interactions between aggregates change 
the functional form of <j>i(z) and Ni(z). In the limit of 
long aggregates 4>i (z) is a linearly decreasing function and 
Ni(z) varies as a square root of height. Although the av- 
erage length of the aggregates can be hardly measured, 
the concentration profile is easily accessible experimen- 
tally and the dependence of <fii(z) can be tested directly. 
The integration of eq. I|13fl gives the expression for the 
global average (Ni) in the range of all heights in the form 
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FIG. 4: Numerical plot of the global average aggregation num- 
ber on the whole range of heights, (Ni) of nomdeal linear 
aggregates, rescaled to its value in the absence of the field, 
(A r i) J 3 = o, as a function of f3 = mgh/kT for different end en- 
ergies: 8 = 15 (thick line, N il0 ^_200), 8 = 25 (thin line, 
JVj, w 2000), 5 = 25 (dash line, N i>0 » 25000). The curve 
was calculated by the integration of eq. 11 -'it for $i = 0.01 
(high concentration). 



m = — n — — (18) 

Similarly to the case of ideal aggregates, (Ni) ~ y ^e 5 
for /3 <C /?* ~ l/JVi,o, while for strong fields it grows up 
with /?, indicating the increasing number of long aggre- 
gates: 



(JVi) = + ^^-/? 2 + 0(/3 4 ) (19) 

It is noteworthy, that the increasing (3 induces not only 
the growth of the aggregates, it also leads to the precip- 
itation of the aggregates to the bottom. At some point, 
the subunit concentration on the top will be close to zero 
and the corresponding average aggregation numbers will 
not be high enough to hold this approximation. To this 
end, the resulting curve of the concentration profile, <f>%(z) 
and the average Ni{z) could be viewed as a combination 
of the above two limits. Namely, Ni(z) ~ y/l — z/z cr 
at the bottom (0 < z < z cr ) and Ni(z) ~ (1 + const 
z)^ 1 on the top (z cr < z < 1). The crossover be- 
tween the regimes, z cr , can be roughly considered as a 
brake up of the approximation of long aggregates (jl7|l . 
Ni(z = z cr ) = 0. That gives^, z cr = v / ^ l n/f3. The con- 
centration profile is presented in Fig. 2 and the average 
length corresponding to that profile depicted in Fig. 3. 
Due to self-assembly, the influence of the gravitational 
field on the concentration profile of living polymers is 



FIG. 5: Three dimensional plot of the distribution function 
of the aggregate lengths at different heights in the isotropic 
phase Pi(N,z) ffity for 8 = 15, (3 = 0.03, = 0.01. All the 
cross-sections with constant heights have the same functional 
form. 



much stronger than that for nonaggregating particles. 
This is an important conclusion since concentration pro- 
files can be measured directly. 

The plot of (Ni) for the whole range of the fields j3 can 
be obtained numerically only. Eqs. fJ-i(z) — const (JSj and 
the normalization condition $j = J Q <f>i(z)dz gives the 
concentration profile as a function of the total subunit 
concentration The profile can be converted into the 
average length as Ni(z) — \J (pi(z)e s and inserted into 
eq. I|13fl . The result is presented in Fig. 4. We plot the 
ratio (Ni) / (Ni)p-o, where (Ni)p = o is its value in the zero 
field. This allows us to compare the relative growth of 
the aggregates with different end energies, 5. As we see, 
the aggregates growth significantly in the gravitational 
field, especially for large S. 

At this point it is important to address the question of 
the magnitude of the gravitational effect for a real sys- 
tem. The typical example of microtubules used in the In- 
troduction gives the value of /3 of the order of 0.001, which 
is significantly lower than the values in Fig. 4. These val- 
ues can correspond to rotating clinostat experiments^ 
performed on microtubules samples, where the system 
rotates about the horizontal axis. The gravitational field 
in these experiments can be up to 300 g. Moreover, the 
influence of the gravity on the global average (Ni) is no- 
ticeable even for normal gravitational conditions in sys- 
tems with low subunit concentrations: taking <f>i = 10~ 5 
and S = 30 we get the average in the absence of the grav- 
ity (A r ,:) /3= o of the order of 10000, while for (3 = 0.001 it 
yields (Ni) ~ 16000, i.e. the normal gravitational field 
induces considerable (60%) growth of the aggregates. 

However, (Ni) cannot show the redistribution of aggre- 
gates of different lengths with height. In order to illus- 
trate the redistribution of aggregates in the gravitational 
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field, it is useful to consider the length distribution func- 
tion Pi(N,z) = J NC l {N,z)dn = 4nNCi(N,z). The eq. 
© gives 



Pi(N,z) = Nexp 



N 



Ni(z) 



(20) 



where we assume iVj(z) in the form (|17|l . The three di- 
mensional plot of this function is given in Fig. 5. It is 
noteworthy, that the envelope curve has the square root 
dependence in accordance with l|17[). The total number 
of aggregates, i.e. the integral of Pi(N, z) on N, is much 
higher at the bottom than on the top. The maximum of 
the curve shifts to higher values of N at the bottom of 
the vessel. 
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IV. NEMATIC ORDERING AND PHASE 
EQUILIBRIUM IN THE GRAVITATIONAL 
FIELD 

Once the subunit concentration exceeds some critical 
value, the nematic phase appears at the bottom of the 
vessel and the system splits into two coexisting phases. 
It is known, that the average length of the aggregates in 
the absence of gravity is different in the nematic and the 
isotropic phasesiSS^S In this chapter we focus on the fol- 
lowing questions: (i) how the gravitational field influence 
the average length of both phases? (ii) is the growth in 
the nematic phase as strong as in the isotropic phase? 
(iii) what would be the distribution function on lengths 
of the aggregates in the isotropic and the nematic phases 
as well as the change of the overall distribution function 
of both phases with increasing gravity? (iv) In particular, 
we will try to answer the question: what is more favorable 
for an aggregate: to grow up in the isotropic phase until 
some critical length and then to change to the nematic 
phase, or to drop out into nematic phase immediately 
in order to decrease its potential energy in the field? In 
any case, it is quite obvious, that the aggregates in the 
isotropic phase will shorten, while the aggregates in the 
nematic phase will grow up with increasing gravity. 

The simplest way to treat the nematic phase is to em- 
ploy the variational method 3 ^*^ that assumes the trial 
orientational distribution function with the adjustable 
variational parameter a, which is found from the mini- 
mization of the free energy In analogy to the case in the 
absence of gravity^^i^ we may write the interaction 
function p n (z) = 4/^/wa(z) and the entropy function 
tr„(z) = ln(a(z)/4) + Na{z)/ (4p), where p is the persis- 
tence length of the rod per the length of one subunit. 3 - 
The last term reflects the flexibility of the rods. The 
only difference with the expression in the absence of the 
fields is that a now depends on the height. Insertion 
of the expressions for a n (z) and p n (z) into J7J) and the 
minimization with respect to P(N, z) along with the con- 
straint J P n (N, z)dN = 4> n (z) gives the expression for the 
distribution function P n (N,z) in the form 



FIG. 6: The global average aggregation number on the whole 
range of heights in the nematic phase, (N n ), eq. I3UL rescaled 
to its value in the absence of the field, (N n )/3 = o, as a function 
of /3 = mgh/kT with different flexibility of the aggregates: 
p = 20 (thick line, N n>0 «_800), p = 50 (thin line, 7V„, w 
600), p = 500 (dash line, N n ,o « 300), for S = 15, $„ = 0.6. 
The overall growth is negligible with respect to the isotropic 
phase, Fig. 4. 



Pn{N,z) 



N 



= exp (N \p- Eoo- j3z- 



P n (N',z) Pn (z)dN' 



a n (z)-S) (21) 



or applying the constraint we can rewrite 



P n (N,z) _ <t> n ( N 

~ivt( 2 ) eXP V N n (z) 



N 



(22) 



where the average aggregation number in the nematic 
phase yields 



N n {z) = S J^<t> n {z)e* 
With this, the free energy of the nematic phase is 



(23) 



F n (z) 

kT 



4>n(z) 



Eoo + (3z + 



a(z) 2 
~ N n {z) 



+ 



7T 4 

-(jy n {z)—== 
4 y/ira(z) 



(24) 



Minimization of the free energy with respect to a{z) 
in the limit of considerably long rods, N n (z) 1, allows 
to relate a(z) and the local concentration as 



a(z) = (2y/np(f> n (z)) 



2/3 



(25) 





FIG. 7: Phase diagram of living polymers in the gravita- 
tional field in the variables: the total subunit concentration 
$ and the parameter associated with the gravitational field 
j3 = mgh/kT for the parameters p — 40, 5 = 10 leading 
to the phase transition even in the absence of gravity. The 
concentrations corresponding to the onset of the phase sep- 
aration diverge and the coexistence region becomes broader. 
Here i denotes the isotropic phase, n is the nematic phase, 
i+n corresponds to the phase coexistence region. 



FIG. 8: Phase diagram of living polymers in the gravitational 
field in the variables: the total subunit concentration $ and 
the end energy of linear aggregates S for p = 100. The dark 
pattern corresponds to the phase separation region in the ab- 
sence of the gravity, /3 = 0, while the light pattern corresponds 
to the gravitational field f3 = 0.2, i denotes the isotropic 
phase, n denotes the nematic phase. The phase separation 
region becomes broader in the gravitational field. 



The chemical potential fi n (z) — dF n (z) / d(j) n (z) is 



kl 



' + ^ n {z)p n {z) (26) 



4p N n {z) 2 

and the osmotic pressure p n (z) = 4>n/ v d [Fn/4>n] /d<j) n (z) 
is 



Pn{z) 

kT 



4>n{z) 



%^) + -^ niZ)Pn{Z \ 



(27) 



The form of these expressions is close to that in the ab- 
sence of gravity^ 

Let us first calculate the average aggregation number, 
N n (z) , in the limit of long aggregates, i.e. when 5 is 
large, so we can neglect the fourth term in l|26l) . Noting 
that ^(f> n (z)p n (z) — a(z)/p, the chemical potential 126|) 



kT 



5 a(z) 

+ l3z+ A 

4 p 



(28) 



Thus, a(z) changes linearly with height, a(z) — a n fi — 
4/5p/3z, where a n fl = a(z = 0). The local concentra- 
tion, is then <f> n (z) = (^ 3 -4/5p 1 / 3 /?z/(2^F) 2 / 3 ) 3/2 

with <p n fi being the concentration at the bottom, and 
correspondingly, the average aggregation number is 



N n {z) 



iV n,0 



5 \8^p 



2/5 



-, 5/4 



p(3z 



(29) 



where N n = N n (z = 0) = ^\/a n (f) n Qe s . Hence, we 
can conclude that N n (z) changes almost linearly with 
height in contrast to the square root behavior in the 
isotropic phase. Since the nematic phase is composed 
mostly of long aggregates, the approximation of long ag- 
gregates works well in a large range of concentrations. 

Assuming, that the vessel is entirely occupied by the 
nematic phase, we can calculate the global average length 
of the aggregates on the whole range of heights 



J N n (z)dz 



(30) 



Although an analytical expression for this case is not 
possible, we can plot the numerical result. Integrating eq. 
(|29|l and fixing the total concentration <!>„ = J Q (f> n {z)dz, 
we get with the help of eqs. and (|2*5|l the depen- 

dence of (N n ) on f3 (Fig. 6). Here again we plot it as a 
rescaled variable, (N n ) / (N n ) p= , in order to compare the 
growth of aggregates with different stiffness, p. As in the 
case of the isotropic phase, the gravitational field leads 
to the growth of the aggregates, though the growth in 
the nematic phase is less pronounced. One can see that 
stiff aggregates (p is large) growth better than flexible 
ones (p is small). By further increasing /3 the top of the 
vessel rarefies until the appearance of the isotropic phase 
coexisting with the nematic phase at the bottom. 

To treat the phase separation region we can employ the 
same method as in ref. |52 which is valid for the coexist- 
ing phases with concentration gradients. Once a system 
is subject to a concentration gradient, there is no concen- 
tration which determines the equilibrium properties of a 
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FIG. 9: Average aggregation number, N(z), as a function of 
height, z, in the phase coexistence region for three different 
parameters associated with the gravitational field, f3 — 0.6 
(thick solid curve), (5 = 1.0 (thin solid curve), j3 = 3.5 (dash 
curve). The jump corresponds to the position of the phase 
boundary, x, the nematic phase is at the bottom, < z < 
x, and the isotropic phase is on the top, x < z < 1. The 
parameters used: p = 40, 5 = 15, $ = 0.255. 



phase. Instead, we will consider the average concentra- 
tion in the nematic phase, which is assumed to be 
at the bottom, z < x, and the average concentration in 
the isotropic phase, $i, at the top z > x. where x is the 
height of the phase boundary: 



- / 4> n {z)dz 



1 



4>i(z)dz 



(31) 
(32) 



The concentrations are related via the conservation of 
mass condition $ = £<£>„ + (1 — x)$i, where <3? is the 
total subunit concentration in both phases. Note, that 
and <&i are not fixed, they are determined from the 
equilibrium conditions. 

The equilibrium corresponds to the equality of osmotic 
pressures and chemical potentials at the boundary, z = x. 
Using the expressions for the isotropic phase (©, i|10|l an d 
the nematic phase, l|26|). i|27|) . we can write 



1 



a(x) 
Ap 



1 



\a(x)(j) n (x)e 5 



^(3)-^= (33) 

<Pn(x) 

^Jja(x)cj) n (x)e s 

E#( X )_4 (34) 
4 y/Tra(x) 



FIG. 10: The position of the boundary between two phases, 
function of /3 — mgh/kT for two different total subunit 
concentrations: $ = 0.255 (thick line) and $ = 0.3 (thin line). 
The parameters used: 5 — 15, p — 40. 



where a(z) is determined from eq. (|25|l . The solution 
of four equations (|31|I - H34|I along with the conservation 
of mass condition, $ = x$„ + (1 — x)&i, gives the equi- 
librium distribution of aggregates between two phases. 
Note, that the last two equations determine the boundary 
concentrations of the coexisting phases, 4>i(x) and <fi n (x). 
The form of the equations coincides with that in the ab- 
sence of gravity^ since the aggregates at the boundary 
have all the same potential energy and the gravitational 
field does not influence their behavior. Chemical poten- 
tials in both phases do not depend on the height. They 
can be regarded as equations for the definition of local 
concentrations 4>i. n {z) as implicit functions of height. To 
get the average concentrations of coexisting phases, $ n 
and <£>i, as well as the boundary position, x, we have 
to integrate the local concentrations. However, it only 
leads to implicit definition of local concentrations. This 
results in a system of integral equations which can be 
hardly solved. See Appendix how to avoid this difficulty. 

In contrast to the phase separation in a homogeneous 
system, the coexisting phases in a gravitational field ex- 
hibit the concentration gradient and so, the average con- 
centrations of the phases, <I>i and <!>„, do not correspond 
to the onset concentrations of the phase transition. This 
is because they both depend on the total concentration, 
<3>. In order to find the values of the total concentration 
signifying the phase separation boundaries, we put in the 
equations x — for the onset of the isotropic phase and 
x = 1 for the onset of the nematic phase. Solving the 
above equations with respect to $, the boundaries of the 
coexistence region can be obtained. Fig. 7 represents the 
typical phase diagram in the variables j3 — The two 
curves diverges, thus, the gravitational field significantly 
broaden the coexistence region. Starting at some point 
in the isotropic phase at zero field and moving vertically 
on the plot, i.e. increasing the field, at some (3 we enter 
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N 

FIG. 11: Length distribution functions in the nematic phase, 
P n (N) (dash curves), and in the isotropic phase, Pi{N) (solid 
curves), for two cases: in the absence of the field, subscript 
j3 = 0, and in the field, subscript (5 = 2.0. The parameters 
used: p = 40, 8 = 15, <E> = 0.255. 



the phase separation region: increasing local concentra- 
tion at the bottom reaches the limiting value and the 
nematic phase appear at the bottom. Increasing gravity 
in the nematic phase leads to the rarefaction on the top, 
what gives rise to the isotropic phase. Thus, the gravity 
facilitates the isotropic to nematic phase transition. 

Since we describe the self-assemble system, which is 
sensitive to the variational of external conditions, in par- 
ticular to the change of the temperature, it is useful to 
plot the phase diagram in the variables 5 — where the 
end energy, S, usually has the inverse temperature depen- 
dence. The result is presented in Fig. 8. Here the inner 
phase diagram (dark pattern) is the coexistence region in 
the absence of the field in accordance with refsi2£*24 As 
we sec, to induce the phase separation, 5 should exceed 
some critical value, for small 6 the aggregates are too 
short and the nematic phase never appears. (Compare 
it with theoretically predicted and experimental phase 
diagrams. The phase diagram depends strongly on 
the stiffness parameter p. In the limit of rigid rods, 
p — ► oo, the concentrations of coexisting phases tend to 
Y-axis, while for flexible aggregates, there is a little gap, 
i.e. stiff aggregates much easily form the nematic phase 
than the flexible ones. The gravitational field severely 
increases the coexistence region (light pattern) in both 
cases. So, for strong enough gravitational fields the left 
border can approach the Y-axis, signifying the absence of 
the pure isotropic phase for high values of S even for very 
low concentrations. This situation should be taken into 
account, when the self-assembled aggregates are subject 
to the ultracentrifugal field. 

Let us focus on the coexistence region in more details. 
The solution of four equations (|3l"1) - (|34fl along with the 
conservation of mass condition for each f3 gives not only 
the average concentrations in the coexisting phases, it 




FIG. 12: The sum of the length distribution functions of the 
nematic and isotropic phases, presented in Fig. 11, P(N) = 
xP n (N) + (1 - x)Pi(N) for (3 = (thick line) and /3 = 2.0 
(thin line). The parameters used: p = 40, 5 = 15, $ = 0.255. 

also gives the position of the phase boundary, x, and the 
concentrations at the boundary, (f>i(x) and 4> n (x). This 
information is enough to obtain the full concentration 
profile composed of two parts, 4> n {z) at < z < x and 
4>i(z) at x < z < 1. This profiles can be converted into 
the average aggregation numbers in both phases: 

7?i(M*)) = \J<t>i(z)e 5 (35) 
NniMz)) = \/(2V^) 2/3 <^ /3 (*) e V4 (36) 

The resulting plot for different (3 is shown in Fig. 9 (Ac- 
tually, we plot Zi >n {(p) and thus, Zj, n (iV) as explained 
in Appendix). The average aggregation number in the 
isotropic phase exhibits the square root behavior on 
height, and it has almost linear dependence in the ne- 
matic phase in accordance to our previous estimates. The 
only exception is the region at the very bottom for high 
gravitational field (3 = 3.5, where N n does not depend 
on z. This is because the concentration at the bottom 
for high gravitational fields is very high and the second 
virial approximation does not work. Thus, after correct- 
ing this approximation, the sharp edge on the dashed 
curve should smooth out. The length of the concentra- 
tion jump at the phase boundary does not depend on the 
field, since the field can only shift the phase boundary 
and it does not influence the equilibrium of the aggre- 
gates lying in the same height. 

As in the case of monodisperse rods of fixed length, 
the gravity can either increase or decrease the volume of 
the nematic phase depending on the total subunit con- 
centration (Fig. 10). If the total concentration is close 
to the isotropic phase boundary and the volume of the 
nematic phase is low, the gravitational field induces the 
increase of the nematic phase, so the volume increases 
until exhausting of the isotropic phase. Since then, the 
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volume gradually decreases indicating the sinking of the 
aggregates in the nematic phase (solid line) . If the total 
concentration is close to the nematic phase boundary, the 
aggregates in the nematic phase can only sink, decreasing 
the volume (thin line). 

Now, let us address the question of the polydispersity of 
the nematic and the isotropic phases in the gravitational 
field. The length distribution functions of aggregates in 
coexisting phases, Pi(N) and P n (N) can be calculated as 



P(N) 

Pn(N) 



N 



1 Uz) 
= 2 exp 



i-xJx N-(z) 
N r <f> n {z) 



N 



x 



° N n (z) 



exp 



Ni(z) 
N 



N n {z) 



dz 



dz (37) 
(38) 



The result is presented in Fig. 11. See Appendix and 
especially eqs. ljA.71) . ljA.81) for the details of the calcu- 
lation. The gravity induces the transition of subunits 
from the isotropic to the nematic phase. So, the area 
of the isotropic distribution function, Pj(JV), decreases 
with increasing gravity, while the area of P n (N) increases 
with p. Nevertheless, the distribution of lengths in the 
nematic phase is much broader in the presence of the 
field, i.e. the gravity increases seriously the polydisper- 
sity of the nematic phase. This is for account of the frac- 
tion of extremely long aggregates induced by the gravity, 
while the fraction of short aggregates remains almost un- 
changed. The isotropic phase looses foremost long, and 
thus, heavy aggregates, but the change in its polydisper- 
sity is not so drastic as in the nematic phase. Hence, 
the precipitation into nematic phase is more favorable 
process than the growth in the isotropic phase. To de- 
scribe the polydispersity of the whole two phase system, 
we can consider the sum of the length distribution func- 
tions, P = xP n (N) + (1 - x)Pi(N). The area of this 
function is proportional to the total number of subunits 
in the system, so its value is the same for any (3. The plot 
of P(N) is shown in the Fig. 12. It brings us to the con- 
clusion that the system as a whole is more polydisperse 
in the gravitational field. Sharply picked function in the 
absence of gravity fade out with increasing gravitational 
field. Note, however, that the maximum of the curve al- 
most does not change. This is the consequence of the 
redistribution of subunits between two phases without 
considerable growth of the aggregates. 



In this paper we focus on the major effects of the grav- 
itational field on the process of self-assembly of linear 
aggregates. It is shown, that the concentration gradi- 
ent induced by the gravitational field leads to the re- 
distribution of aggregates of different lengths (and thus, 
masses) on different heights: the local concentration of 
the self-assembled rods interacting via the second virial 
coefficient decreases linearly with height in contrast to 
barometric distribution of aggregates with fixed length. 
This prediction can be directly checked experimentally 
in the rotating clinostat experiments. In turn, the aver- 
age length of the aggregates has a square root depen- 
dence in contrast to hyperbolic dependence found for 
ideal aggregates^ We also show, that except the redis- 
tribution of aggregates with height, the gravity induces 
the overall growth of the aggregates. The stiffness of the 
aggregates enhance the growth under gravity. Strong 
enough gravitational field can provoke the isotropic to 
nematic phase transition, which is reflected in the sig- 
nificant broadness of the phase coexistence region. The 
transition manifests itself in the precipitation of the ag- 
gregates from the isotropic to the nematic phase regard- 
less their length. Furthermore, the length distribution of 
the aggregates is found to be very sensitive to the gravita- 
tional field: the gravity induces considerable polydisper- 
sity in the nematic phase. Taken together, the effect of 
gravity can lead to a series of experimentally relevant ob- 
servations which are interesting from fundamental point 
of view. 



V. CONCLUSION 



Although the gravitational field is often neglected in 
the description of physicochemical processes, it may have 
noticeable influence on such biological heavy objects as 
microtubules or actin filaments, essential components of 
living cells. The equilibrium properties of such objects is 
influenced even by the earth gravitational acceleration. 
The fields in an ultracentrifuge can multiple the effect of 
gravity. 
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APPENDIX: ' 

Equations u in = const (jSJ), i|26|) do not allow to obtain 
the analytical expression of local concentrations 4>i_ n {z) 
as a function of z except for two asymptotes eqs. (II 1 1) . 
(|16fl . To this end, the total concentration, $ = J <fi(z)dz, 
and the average lengths, Ni }Tl {z), are, in fact, the inte- 
gral equations for the functions <f>i^ n {z), the functional 
forms of which are determined from the chemical poten- 
tials. In general, it is rather difficult to solve such sys- 
tem of equations. However, we can employ the structure 
of equations ©, (|26[> . which allows to get analytically 
the inverse functions, Zi tn ((f>). Namely, for the isotropic 
phase 



where the value of the local concentration is </>„ (0) at the 
bottom and <fii(l) on the top. 

The derivatives in the integrals are 



dzi{4>) 



dz n ((j>) 



2/3 
5 



(A.5) 



, ^ (W V2 + 

3 (2V^P) 7 V 1/6 /? V 1 



Thus, the average concentrations can be integrated an- 
alytically. 

The global length average of the aggregates in both 
phases, (|13[) and ( I3U|) . can be written as integrals on <f>: 



Zi{<t>) =x - - 



1 



\f<t>i{x)t 



4>e 5 * 



(A.l) 

Here we assume the isotropic phase is situated above 
the phase boundary, x < Zi < 1, and the boundary con- 
centration is <pi(x). In case of the homogeneous phase, 
one can put in this expression x — 0. Analogously, the 
inverse function for the nematic phase yields in the form 



(JVi) = ^e s ' 2 II ^1 



dej) 



(20Fp) 1/3 / V „ 



,(0) 



(A.7) 



# (A.8) 



The length distribution functions l|37|) . (|38() yield in the 
form 



Zn{4>) = X + 

4 



2 4 / 3 7r 1 /6pl/3 e <5/2 /3 

(\/4>^-i/^\x)) 

2 4 / 3 7r 1 /6pl/3 e 5/2 /3 



+ Pi(N) = 



(A.2) P n (N) = 



The two expressions should be accompanied by the 
boundary conditions, Zi((p = </>i(l)) = 1 and z n {4> = 
^n(O)) = 0. 

With this, the average concentrations in the phases are 



1 Bz.^) 



l-X 



<t>i{*) 



1 f*"W dZn(& 



(A.3) 
(A.4) 



AT 
l-x" 



,-5 



0,(1) 



exp 



\ dz^ 



AjB) 



N Ae- S 




<t> n (x) 



*„(0) 

2JV 



dz n {cj>) 



4.10) 



where we used (|35ll , (|36|) for the averages in the isotropic 
and the nematic phases. 

Solution of algebraic eqs. (j3HJ), (|A.3|I . ( IA.4|I along 
with the conservation mass condition, $ = x$ n + (1 — 
gives the phase diagram as well as the equilibrium 
average concentrations in the coexisting phases. 
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